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Abstract 



This Letter deals with topological solitons in an 0(3) sigma model in three space dimensions 
(with a Skyrme term to stabilize their size). The solitons are classified topologically by their Hopf 
number N. The N = 2 sector is studied; in particular, for two solitons far apart, there are three 
"attractive channels" . Viewing the solitons as dipole pairs enables one to predict the force between 
them. Relaxing in the attractive channels leads to various static 2-soliton solutions. 



1 Introduction 



The 0(3) nonlinear sigma model in 3+1 dimensions involves a unit vector field <j> = {4 ,1 >4 ,2 i4> 3 ) 



which is a function of the space-time coordinates x 11 = (t, x, y, z). Since <j) takes values on the unit 
sphere S 2 and the homotopy group ^(S" 2 ) is non-trivial, the system admits topological solitons 
(textures) with non-zero Hopf number N G ^(S* 2 ). The model arises naturally in condensed-matter 
physics and in cosmology, and its dynamics in these contexts is governed by the Lagrangian density 
(dfjcf)) 2 . With such a Lagrangian, there are no stable soliton solutions: the usual scaling argument 
shows that textures are unstable to shrinking in size. (In practice, the decay of these textures is 
more complicated, and can, for example, involve decay into monopole-antimonopole pairs.) 

Another context in which the nonlinear sigma model arises is as an approximation to a gauge 
theory; ie one gets an effective action for the gauge field in terms of the scalar fields 4> a . In this case, 
there are higher-order terms in the Lagrangian, which can have the effect of stabilizing solitons. 
This idea is familiar in the Skyrme model, where the target space is a Lie group; it also occurs for 
the case relevant here, where the target space is S 2 |lj , [|] , || . 

The simplest such modification involves adding a fourth-order (Skyrme-like) term to the La- 
grangian. This leads to a sigma model which admits stable, static, localized solitons — they 
resemble closed strings, which can be linked or knotted. The system has been written about at 
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least since 1975 Q, [||; but recently interest in it has increased, stimulated by numerical work ||, 
§, §, |, @, @; see also @, @, Q. 

In this Letter, we shall deal only with static configurations, so <j> is a function of the spatial 
coordinates x 3 = (x,y,z). The boundary condition is (j> a — ► (0,0,1) as r — > oo, where r 2 = 
x 2 + y 2 + z 2 . So we may think of as a smooth function from S 3 to S 2 ; and hence it defines a 
Hopf number TV (an integer). This N may be thought of as a linking number: the inverse images 
of two generic points on the target space, for example (0,0, 1) and (0,0, —1), are curves in space, 
and the first curve links N times around the other. 

The energy of a static field <f){x^) is taken to be 

E ^3^J WifWiF) + FjkFjk] d 3 x, (1) 

where Fjk — e a bc4' a (dj(j) b )(dk4> c )/2. The ratio of the coefficients of the two terms in (|l|) sets the 
length scale — in this case, one expects the solitons to have a size of order unity (note that other 
authors use slightly different coefficients). The factor of 1/327T 2 is justified in [fl2|| : there is a lower 
bound on the energy which is proportional to N 3 / 4 , and if space is allowed to be a three-sphere, 
then there is an N = 1 solution with E = 1. So one expects, with the normalization (|l]), to have 
the lower bound E > N 3 / 4 . 



2 The One-Soliton 

The minimum-energy configuration in the N = 1 sector is an axially-symmetric, ring-like structure. 
It was studied numerically in [[| (with no quantitative results) , in [Q (which gave an energy value 
of E = 1.25, although without any statement of numerical errors), and in 0, [[To] (where the field 
was placed in a finite- volume box, so its energy was not evaluated accurately). For the results 
described in this letter, a numerical scheme has been set up which 

• includes the whole of space R 3 (by making coordinate transformations that bring spatial 
infinity in to a finite range); and 

• using a lattice expression for the energy in which the truncation error is of order /i 4 , where h 
is the lattice spacing. 

Using this shows that the energy of the one-soliton is E = 1.22 (accurate to the two decimal places). 

Let us choose the axis of symmetry to be the z-axis, and the soliton to be concentrated in the 
xy-plane. In terms of the complex field 

(the stereographic projection of </>), the 1-soliton solution is closely approximated by the expression 

x + iy 
z - if{r) 
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where f(r) is a cubic polynomial. We may minimize the energy of the configuration with respect 
to the coefficients of /: this gives an energy E — 1.23 (less than 1% above the true minimum), for 

f(r) = 0.453(r - 0.878)(r 2 + 0.705r + 1.415). (4) 

Note that W — > as r — > oo (the boundary condition); that W = on the z-axis; and that 14 7 = oo 
on a ring (of radius 0.878) in the rry-plane. The ring where W = oo links once around the "ring" 
(the z-axis plus a point at infinity) where W — 0: hence the linking number N equals 1. The field 
looks like a ring (or possibly a disc, depending on what one plots) in the xy-plane. 

To leading order as r — > oo, 1 and 4> 2 have to be solutions of the Laplace equation (cf) 1 and cf> 2 
are the analogues of the massless pion felds in the Skyrme model). From (Q) one sees that 

cP 1 + i<f> 2 « 2W « for larger. (5) 

So (ft 1 and <j) 2 resemble, asymptotically, a pair (P, Q) of dipoles, orthogonal to each other and to 
the axis of symmetry. 

It is useful to note the effect on the 1-soliton field of rotations by 7r about each of the three 
coordinate axes. These, together with the identity, form the dihedral group £>2- Let ±7 and ±C 
denote the maps 

± I : W i-> ±W, ±C : W i-> ±W. (6) 

Then it is clear from (||) that the four elements of Di induce the four maps {/, —I, C, —C} on W. 

The single soliton depends on six parameters: three for location in space, two for the direction 
of the z-axis, and one for a phase (the phase of W). The phase is unobservable, in the sense that 
it does not appear in the energy density, and can be removed by a rotation of the target space S 2 ; 
but the energy of a two-soliton system depends on the relative phase of the two solitons, as we shall 
see in the next section. 

3 Two Solitons Far Apart 

Suppose we have two solitons, located far apart. Let (P + ,Q + ) denote the dipole pair of one of 
them, (P_,Q_) the dipole pair of the other, and R the separation vector between them. There 
will, in general, be a force between the solitons, which depends (to leading order) on the distance 
R between them, and on their mutual orientation. One can predict this force by considering the 
forces between the dipoles. Since the fields are space-time scalars, like charges attract; so the force 
between two dipoles is maximally attractive if they are parallel, and maximally repulsive if they 
are anti-parallel. There are three obvious "attractive channels" (mutual orientations for which the 
two solitons attract), which will be referred to as channels A, B and C. We now discuss each of 
these. 

Channel A. The only axisymmetric configuration involving two separated solitons is one where 
each dipole pair is orthogonal to the separation vector: in fact, where P+ x Q + , _P_ x Q- and 
R are all parallel. The configuration is illustrated in Fig 1. The two circles are where W = oo, 
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Figure 1: Two solitons in channel A, and the corresponding static 2-soliton solution. 



while the line linking them is where W = 0. The arrows on the curves serve partly to distinguish 
solitons from anti-solitons: the convention here is that solitons obey the right-hand rule, whereas 
anti-solitons would obey the left-hand rule. 

Let 9 denote the angle between P+ and P- (ie the pair (P+, Q+) is rotated by 9 about the line 
joining the two solitons). Let E\ denote the energy of a single soliton, and E 2 (R, 9) the energy of 
the two-soliton system, as a function of the separation R and the relative phase 9. Considering the 
potential energy of the interacting dipoles suggests that 

E 2 {R,9) = 2Ei - 2kR- 3 cos9 (7) 

for some constant k. Clearly E2(R,9) is minimized, for a given R, when 9 = (ie when the two 
solitons are in phase): this is channel A. The formula (^) was tested numerically, by computing 
the energy of the configurations obtained by combining translated and rotated versions of the 
approximate one-soliton (||) . The combination anstaz was simply that of addition (W = W+ + W- ) , 
which is a plausible approximation for large R (bearing in mind that the VF-field tends to zero away 
from each soliton). For R in the range 6 < R < 16, the form (Q) is indeed found to hold, with k w 1. 
(In view of the crudity of the "sum" ansatz, the accuracy is not claimed to be better than 10% or 
so; but the R~ 3 cos9 dependence is very clear.) The behaviour under the discrete symmetries D 2 
is the same as for the 1-soliton, namely {/, —I, C, — C}. 

Channel B. This channel is one in which both dipole pairs are co-planar with R, with P + x Q + 
and P- x Q- being parallel (and orthogonal to R). This configuration is depicted in Fig 2(a). In 
this case, the effect of the discrete symmetries D2 on the configuration is {/, /, C, C}. Consideration 
of the forces between the dipoles suggests that the energy behaves like 

E 2 (R, 6) = 2Et + kR- 3 cos(6» + - 6 J), (8) 

where 9± is the angle that P± makes with R. The expression (||) has a minimum when 9 + — 9- = it, 
and this is attractive channel B. 
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Figure 2: Two solitons in attractive channel B coalesce to form a static 2-soliton solution. The pictures 
(a), (b) and (c) are a time-sequence. 




Channel C. Here the dipole pairs are again co-planar with R, but now P + x Q + and P_ x Q_ 
are anti-parallel. This is depicted in Fig 3(a). In this case, one expects that 

E 2 {R,6) = 2Ei + 3/ciT 3 cos((9+ + 6L), (9) 

where, as before, 8± is the angle that P± makes with R. So the attractive force is maximal 
when 6 + + 6- = 7r. Fhe dependence was confirmed numerically, as before, with This 
'maximally-attractive' channel is referred to as channel C . Fhe effect of the discrete symmetries 
D 2 is {I, I, C, C}, as for channel B. 

4 Relaxing in Channel A 

In this section, we see what happens when we begin with two solitons far apart (in the first of the 
attractive channels described above), and minimize the energy. Fhis was done numerically, using 
a conjugate-gradient procedure. 

Suppose, then, that we start in the (axisymmetric) channel A, and minimize energy without 
breaking the axial symmetry. Fhen the two solitons approach each other along the line joining them, 
and the minimum is reached when they are a nonzero distance apart. Fhe resulting configuration 
therefore is a static solution of the field equation; it has energy E — 2.26 (this is to be compared 
with 2E\ = 2.45), and it resembles two rings around the z-axis, separated by a distance R = 1.3. In 
other words, W = consists of the 2-axis plus infinity, as for the 1-soliton, while W — 00 consists 
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Figure 3: Two solitons in attractive channel C. 




of two disjoint rings around it; so the linking number N does indeed equal 2. The picture is as in 
Fig 1. 

There is an approximate configuration analogous to (^), namely 

w _ 2{x + iy)[z{l + I3r) - ih{r)] 

[z + (l + ar + )-if(r + )][z-(l + ar-)-if(r-)Y [ 1 

where a and (3 are parameters, / and h are cubic polynomials, z± — z±0.65 and r± — x 2 +y 2 + z±. 
The two- ring structure is evident from ( |Ic| ) . Minimizing the energy of (|l(]) with respect to the ten 
parameters (a, (3, and the coefficients of / and h) gives an energy E = 2.29, it ie 1.3% above that 
of the solution. The corresponding configuration <j) is very close to the actual solution. 

While this is a solution, it is not the global minimum of the energy in the N — 2 sector; in 
particular, channel B produces a solution with lower energy. So the question arises as to whether 
the channel A minimum is stable to (non-axisymmetric) perturbations {ie whether it is a local 
minimum of the energy, as opposed to a saddle-point). The linking behaviour of the channel B 
minimum is that of a single ring around a double axis (as we shall see in the next section), as 
opposed to a double ring around a single axis; there is a continuous path in configuration space 
from the one configuration to the other, but the contortions involved in this suggest that there is 
an energy barrier (in other words, that the channel A solution is a local minimum). Numerical 
experiments, involving random perturbations of this solution, provide strong support for this; but 
more study is needed. 



6 



5 Relaxing in Channel B 



Next, we start in channel B and once again flow down the energy gradient. As depicted in Fig 2, 
the two rings (where W = oo) merge into one, and then the two lines where W = merge as well. 
We end up with a solution which has been described previously @ , [(§ , Jl(| , and which is believed 
to be the global minimum in the N — 2 sector. It is axially-symmetric, and resembles a single ring; 
but this time the ring winds around a double copy of the z-axis, and hence it has a linking number 
of N = 2. The energy of the solution is E = 2.00, which agrees with the figure given in 0. 

As before, we can write down an explicit configuration which is very close to the solution. One 
such expression is 

W=^%, (ID 
azr — ij(r) 

where a is a constant and f(r) is a quintic polynomial. Minimizing the energy with respect to the 
six coefficients contained in ( |ll| ) gives E — 2.03 [ie 1.5% above the true minimum), for 

a = 1.55, f(r) = 0.23(r - 1.27) (r + 0.44) (r + 0.16)(r 2 - 2.15r + 5.09). (12) 

Since / has only one positive root, W — oo is a ring (of radius 1.27) in the xy-plane; whereas 
W = is the z-axis, with multiplicity two. The components of <p derived from (|ll| ) are very close 
to those of the actual solution. 



6 Relaxing in Channel C 

If one begins with the configuration depicted in Fig 3(a) and moves in the direction of the energy 
gradient, the two solitons approach each other. If the two W — oo loops touch, one has a figure-eight 
curve, with the W = lines linking through it in opposite directions: Fig 3(b). This configuration is 
certainly not stable: preliminary numerical work indicates that the two 'halves' of the configuration 
rotate by 7r/2 (in opposite directions) about the axis joining them. So the figure-eight untwists to 
become a simple loop, and the two W = curves end up pointing in the same direction, exactly as 
in Fig 2(b) and (c). Hence the minimum in channel C is the same as that in channel B. Between 
this mimimum and the channel- A one, there should be saddle-point solutions; but what these look 
like is not yet clear. 



7 Concluding Remarks 

There has already been some study of two-soliton dynamics, using a "direct" numerical approach 
(see, for example, 0); this is computationally very intensive. The results reported in this Letter 
could be viewed as the first step towards a somewhat different approach, namely that of constructing 
a collective-coordinate manifold for the two-soliton system. The analogous structure for the Skyrme 
model has been investigated in some detail [Est , [fliif ; in particular, it has the advantage that one can 
introduce quantum corrections by quantizing the dynamics on the collective-coordinate manifold 
p%. Since each Hopf soliton depends on six parameters, the two-soliton manifold M 2 should have 
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dimension (at least) twelve; each point of M2 corresponds to a relevant N = 2 configuration, and 
the expressions (||) and (|ff|) are examples of such configurations. 

But clearly much more work remains to be done towards understanding the energy functional 
on the N = 2 configuration space. The suggestion of this Letter is that the global minimum (which 
is, of course, degenerate: it depends on six moduli) is as in Fig 2(c); there is a local minimum as 
in Fig 1; and between the two are saddle-point solutions which may be related to the figure-eight 
configuration Fig 3(b). 
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